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Abstract. We investigate the characteristics and the time evolution of the cosmic web from 
redshift, 2 = 2, to present time, within the framework of the NEXUS+ algorithm. This ne¬ 
cessitates the introduction of new analysis tools optimally suited to describe the very intricate 
and hierarchical pattern that is the cosmic web. In particular, we characterize filaments (walls) 
in terms of their linear (surface) mass density. This is very good in capturing the evolution of 
these structures. At early times the cosmos is dominated by tenuous filaments and sheets, which, 
during subsequent evolution, merge together, such that the present day web is dominated by 
fewer, but much more massive, structures. We also show that voids are more naturally described 
in terms of their boundaries and not their centres. We illustrate this for void density profiles, 
which, when expressed as a function of the distance from void boundary, show a universal pro¬ 
file in good qualitative agreement with the theoretical shell-crossing framework of expanding 
underdense regions. 
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1. Introduction 

On megaparsec scales the matter distribution of the Uni verse is n ot uniform , but it 
forms an intricate pattern which is known as the Cosmic Web (IBond. Kofman fe Pogosvan 
1996t h This is the most salient feature of the anisotropic gravitational collapse of matter, 


the motor behind the formation of structure in the Universe. Identifying and character¬ 
ising the cosmic web network, in both numerical simulations and observations, is very 
challenging due to the overwhelming complexity of the individual structures, their con¬ 
nectivity and their intrinsic multiscale nature. It is even more difficult to follow the time 



Figure 1 . The complexity and multiscale character of the cosmic web as identified by 
NEXUS+. It shows the density field (left), the filaments (centre) and the walls (right) in a 
100 x 100 x 10 (/i -1 Mpc) 3 slice trough the Millennium-II Simulation. 
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evolution of the cosmic web, since the dominant components and scales change rapidly 
with redsliift. This necessitates the use of scale- and user-free methods that naturally 
adapt to the complex geometry of the web and that extract the maximum information 
available among the components of this network. 


2. Simulation and methods 

We follow the evolution of the cosm ic web using the high resolution Millennium an d 
Millennium-II dark matter simulations ( Sprineel et all2005t Bovlan-Kolchin et al.ll2009h , 
which describe the formation of structure in a periodic box of length 500 a nd 100 fe~ 1 Mpc, 
respectively. We use the Delaunay Tessellation Field Estimator (ISchaan fe van de Wevgaert 


200Ci van de Wevgaert. & Schaan 2009HCautun fc van de Wevgaerf} 2011 1 to obtain con- 


tinuous density and velocity fie lds. These are used as input for the NEXUS+ algorithm 
ijCautun. van de Wevgaert & Jones 201 3). which identifies the cosmic web components. 
NEXUS+ uses the mass distribution, as traced by the density field, to identify clusters, 
filaments, walls and voids. The dominant morphological signal is extracted from a 4- 
dimensional scale-space representation of the matter distribution. The fourth dimension 
is constructed by filtering the logarithm of the density with a Gaussian kernel, for a range 
of smoothing scales. The environments are classified on the basis of the dominant local 
morphological signature of the filtered density field. 

The outcome of applying NEXUS+ to the density fiel d is illu strated in Fig. [1] Based 
on the Multiscale Morphology Filter (MMF) method of Aragon-Calvo et al. (12007b} ). we 
developed NEXUS+, and its sister method, NEXUS, to use a multitude of tracer fields 
for classifying the cosmic web environments. It makes use of the density, tidal, velocity 
shear and velocity divergence fields. In addition, NEXUS and NEXUS+ have the unique 
feature of employing natural and self-consistent criteria for identifying the cosmic web, an 
improvement with respect to the less clear percolation threshold criteria used in MMF. 
This leads to approaches optimally suited for identifying the morphological environments 
since these methods are multiscale, parameter-free and designed to fully account for the 
anisotropic nature of gravitational collapse. 

For this study, NEXUS+ has two main advantages. Firstly, it determines in a self- 
consistent way all the different morphological components. And, secondly, the multiscale 
character of the method makes it ideal not only for identifying prominent as well as 
tenuous environments, but also for studying the time evolution of the cosmic web without 
having to choose a user-defined scale. 


3. The evolution of the cosmic web 

Fig. [2] shows the evolution of filamentary environments starting with a redshift of 
z = 1.9 down to the present time. At early times, the filaments form a complex network 
that pervades most of the cosmic volume, with the exception of the most underdense 
regions. While the network has a few thick structures, it is dominated by small scale 
filaments. These thin filaments seem to be packed much more tightly close to prominent 
structures, suggesting that overdense regions have a higher richness of filaments. By 
z = 1, we find that most of the tenuous structures have disappeared and that we can 
more easily see the pronounced filaments. Going forward in time, to z = 0.5 and 0, 
we find that the evolution of the cosmic web significantly slows down, with only minor 
changes after z = 0.5. Though not shown, a simila r ev olution can be seen for the wall 
network too. For a more in depth analysis see lCautun et al.1 ( 2014 ). 
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Figure 2. The evolution of the filamentary network as identified by NEXUS+ in a 
100 x 100 x 10 (/i -1 Mpc) 3 slice trough the Millennium-II Simulation. 


3.1. The mass distribution along filaments and walls 

One of the most widely employed methods to study the cosmic web properties involve 
the use of global properties, like mass and vo l ume filling fraction of each component (e.g. 
Araeon-Calvo et al. 2007bt lHahn et al' 2007 : Ihbrero- Romero et al. 2009ll . When applied 


to our study, such an approach leads to the results presented in top row of Fig. [3] For 
example, it shows that the volume occupied by filaments decreases since high redshift. 
This is probably due to the merging of the thin and tenuous filaments with the more 
prominent structures. However, such a simple analysis cannot characterise the complex 
evolution seen in Fig. [2] For example, it cannot tell which components, tenuous or promi¬ 
nent, contained the most mass and how this mass distribution evolves in time. To do so, 
one needs a more complex analysis framework that takes into account the geometry of 
the various cosmic components.__ 


Such a framework has been introduced inlCautun et alJ (2014. see alsolAragon-Calvo. van 


2010|) and takes advantage that, to a first approximation, filaments and walls can be seen 


as 1-dinrensional lines and 2-dimensional surfaces, respectively. These represent the spine 
of filaments and the central plane^of sheets; and can b e com puted using the techniques 
described in Aragon-Calvo. van de Weveaert fe Jones! ( 2010l) and Cautun et al. ( 2014ll . 
This reduces the complex filamentary network to a simpler distribution of interconnected 
curves, in which all the mass and haloes in the filament is compressed to its spine. Fol- 
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filament linear density C ( h 1 M @ / [h 1 Mpc] ) wall surface density o ( h 1 M @ / [If 1 Mpc ] 2 ) 


Figure 3. Top row: the time evolution of the mass and volume fraction in each cosmic web 
environment. Bottom row: the distribution of linear mass density of filaments (left) and the 
surface mass density of walls (right) at different redshifts. 


lowing this, one can move along the resulting curves and compute local qua ntities, like 


the mean linear density and the mean diameter of filaments (for details see Cautun et al 

l201ll) . 


The bottom-left panel of Fig. [3] shows the results of such an analysis. It gives the length 
of filaments, per unit volume, that have a given linear mass density, when measured within 
a window of 2 h -1 Mpc. The figure shows that there are few very low or very high mass 
filaments and that most of the length of the filamentary network is given by segments 
with linear densities of ^10 12 Mq/Mpc. This distribution evolves in time, to show that 
segments with high mass become more common at late times, while at the same time 
there are fewer low mass segments. Even more telling, is the shift to the right in the peak 
of the distribution, showing that time evolution leads to an increase in the mass of the 
typical filament segment. 

In contrast to filaments, the bottom-right panel of Fig. [3] shows that the typical sheet 
regions become less massive at present time. The decrease in wall surface density is seen 
as the shift in the peak of the distribution towards lower er wa ii values at later times. It 
shows that the decrease in the mass fraction of walls seen in the top-left panel takes 
place via two processes. First, as we just argued, typical sheet stretches become less 
massive. And secondly, the extent of the wall network reduces at later times, as seen in 
the decreasing peak values of the er wa n distribution. 
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Figure 4. The fraction of haloes in each cosmic web environment as a function of halo mass. 
We shows results for redshift, z = 2 (left), and for the present time (right). 



Figure 5. The haloes populating a few typical cosmic web filaments (left) and walls (right). 
The colour and size of the haloes indicates their mass (see legend on the right) and the black 
points show haloes found in node environments. 


3.2. The halo distribution among the web environments 
The distribution of haloes across web environments plays a key role. First, the mor¬ 
phology of the cosmic web determines the preferential directions of accretion and thus 
influ ences the shape and angular momentum of h aloes, inducing large scale alignments 
(e.g. I Aragon-Calvo et al. 2007a : Hahn et ahl 2007l l. And secondly, to identify and study 
the cosmic web in the observational reality one makes use of galaxies, which are hosted 
in haloes. So knowing how haloes populate the morphological components is crucial to 
understand the web environments seen in observations. To this end, we show in Fig. [4] 
the distribution of haloes in the cosmic web at both z = 2 and 0. The figure shows a clear 
segregation of haloes across environments, with the most massive ones living in nodes 
and prominent filaments. Walls typically host ~40 12 and lower mass objects, 

while void regions are populated with even lower mass haloes. 

More interestingly, the halo population varies not only with environment, but also 
between structures with similar morphological features. This is shown in the left panel 
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Figure 6. The mass transport across the web components. It shows dark matter particles 
identified at a = 2 as part of filaments (left), walls (centre) and voids (right). These particles 
are coloured according to their web classification at z = 0: node (light-grey), filament (grey), 
wall (dark-grey) and void (black). 


of Fig. [5j where we illustrate the haloes found in several filamentary branches. From top 
to bottom, we show progressively less prominent filamentary branches, with the bottom 
most ones corresponding to filaments found in voids. The panel shows a clear trend 
of the halo distribution with filament properties. Thicker filaments, which are typically 
outstretched between cluster pairs, are populated with more massive haloes which are 
also more tightly packed together. In contrast, haloes in tenuous filaments are typically 
low mass, similar to the ones in walls, and are widely spaced apart. 

Fig. [5] explains why prominent filaments are easy to find observationally: they host 
many bright galaxies. In contrast, tenuous structures are mostly inhabited by a dilute 
population of low luminosity galaxies. Detecting these structures presents many obser¬ 
vational challenges since such tenuous objects are hardly conspicuous in the spatial dis- 
tribution of galaxie s. The configuration of three aligned galaxies in side a void found b y 
Bcvku et al.! ( 20131) is probably an example of such a thin filament ( Rieder et al. 2013 ). 


3.3. The mass transport across the cosmic web 
By comparing the cosmic web identified at different redshifts, we can trace the path 
taken by the anisotropic collapse of matter and study the transport of mass among 
the different web environments. According to the gravitational instability theory, the 
matter distribution follows a well defined path, with mass flowing from voids into sheets, 
from sheets into filaments and only in the last step into the cosmic nodes. We illustrate 
this with the help of Fig. [G] which shows the distribution of dark matter particles in 
a thin 2 /i^ 1 Mpc slice at redshift z = 2. Each panel gives the particles identified as 
part of filaments, walls and voids at z = 2, with the particles coloured according to the 
environments they are found in at the present time. 

Fig. [G] shows some of the most important characteristics of the mass transport across 
the cosmic web . T hese conclus ions are supported by the more in-depth and quantitative 
analysis of lCautun et al.1 ( 2014 ). Among others, it shows that nodes form at the inter¬ 
sections of the filamentary network and that most of the mass in them has been accreted 
from regions that correspond to z = 2 filaments. For filaments, most of the mass found 
in these objects at z = 2 is also found in filaments at present time. In contrast, walls 
loose more than half of their z = 2 mass to mostly filaments, with the rest remaining in 
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Figure 7. The density profile, 1 + <5, of voids for the full population and for subsamples split 
according to the volume of the void, from small to large. The left panel show the spherical 
profile of voids. The right panel shows the void profile as a function of the distance from the 
void boundary. We use the convention that negative d values correspond to the inside of voids, 
while positive ones correspond to the profile outside the void. 


walls. And finally, voids also loose around half of their z = 2 mass, which flows into walls 
and filaments. In addition, Fig. [S] illustrates some of the limitations of the analysis. A 
small fraction of z = 2 filaments and walls are identified as present day walls and voids, 
respectively. This is restricted to minor filaments and walls and it is due to the difficulty 
of identifying such tenuous structures, especially across multiple time steps. 


4. A natural void profile 

The filaments and walls are not the only ones that have complex shapes and mor¬ 
phologies. Voids have it too. The simple picture of an expanding underdensity in a 
unifo r m background suggests that vo ids become more spherical as they evolve (llcke 
1198-4 ISheth fc van de Wevgaert 2004 1. But in reality, voids are not isolated and this 
simple picture does not hold. There are two major factors that affect the evolution of 
voids. Firstly, since voids are nearly empty the force field that dictates their growth is 
dominated by external forces. This external tidal field determines the anisotropic ex¬ 
pansion of voids and manifest itself as large scale corr elations, over distances larger than 
30 h~ 1 Mpc, between the shapes of neighbouring voids ( Platen, van de Wevgaert fe Jones 
200 8). Secondly, as voids expand, they encounter neighbouring voids resulting in a pack¬ 


ing. problem. T hese effects l ead to v oids that have highly aspherical and complex shapes 
( Platen, van de Wevgaert fe Jonesl 2007112008 1. 

Up to now, voids were characterised in terms of spherical profiles with respect to the 
void’s center, typically the barycentre. Such a methodology results in the void profiles 
shown in the left panel of Fig. [71 These pr ofiles have been rescaled following the prescrip¬ 
tion of Hamaus. Sutter fc Wandeltl ( 2014 ). where i? e ft is the effective radius of the void 
corresponding to a spherical void with the same volume as the real object. This results 
in void profiles that are not fully universal and where one needs to follow a complex pro¬ 
cedure to rescale void s of different sizes (for details see lHamaus. Sutter fe Wandeltll2014 : 
Nadathur et al.ll2014 ). More worryingly, the resulting profile does not show the large den¬ 
sity ca ustic at their boundary res ulting from the shell-crossing of expanding underdense 
shells ( Sheth fe van de Wevgaertll2004l) . So what is the reason for the mismatch? 

The discrepancy b etween sp herical void profiles and the shell-crossing predictions of 
Sheth & van de Weygaert ( 20041 ) arise from the fact that voids are very far from having 


a spherical shape. Thus, using spherically averaged profiles does not lead to an accu- 
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rate description of void structure. In fact, it is easier to determine the void bound¬ 
ary, where most of the mass and galaxies reside, than the void center, which is devoid 
of tracers. In fact, t his ver y fact is the co rnerst o ne o f Watershed-based void finders 
I Platen, van de Wevgaert fc .Tonesl 2007t Nevrinckl 2008 ). Thus, it is more natural to de¬ 
scribe voids with respect to their boundary than with respect to their center. 

It suggests that void profiles, and in general void properties, should also be computed 
with respect to the void boundary. For this, we define the boundary distance field d ? ; (x) 
as the minimum distance between point x and the boundary of void i. In addition, to 
distinguish between points inside and outside the void, we take the convention that di(x) 
is negative if x is inside void i and positive otherwise. The resulting density profiles, as 
a function of the void boundary distance, are shown in the right panel of Fig. |T| First, 
we find the sharp increase in the density profile at 0, corresponding to the caustic 
resulting from shell-crossing. And secondly, voids of different size show a much more 
similar and universal profile. This method gives a better and more natural description of 
not only density profiles, but also of void velocity profiles (Cautun et al. 2015, in prep.). 
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